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The Laplacian gauge is a nonperturbative gauge fixing that reduces to Landau gauge in the 
asymptotic limit. Like Landau gauge, it respects Lorentz invariance, but it is free of Gribov copies; 
the gauge fixing is unambiguous. In this paper we study the infrared behavior of the lattice gluon 
propagator in Laplacian gauge by using a variety of lattices with spacings from a — 0.125 to 0.35 fm, 
to explore finite volume and discretization effects. Three different implementations of the Laplacian 
gauge are defined and compared. The Laplacian gauge propagator has already been claimed to be 
. insensitive to finite volume effects and this is tested on lattices with large volumes. 
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I. INTRODUCTION 

(N 

. The lattice provides a useful tool for studying the gluon propagator because it is a first principles treatment that 
can, in principle, access any momentum window. There is tremendous interest in the infrared behavior of the gluon 
propagator as a probe into the mechanism of confinement jjj and lattice studies focusing on its ultraviolet behavior 
have been used to calculate the running QCD coupling ||. Such studies can also inform model hadron calculations ||. 
Although there has recently been interest in Coulomb gauge and generic covariant gauges ||, the usual gauge for 
these studies has been Landau gauge, because it is a (lattice) Lorentz covariant gauge that is easy to implement 
on the lattice, and its popularity means that results from the lattice can be easily compared to studies that use 
different methods. Finite volume effects and discretization errors have been extensively explored in lattice Landau 
gauge ||. Unfortunately, lattice Landau gauge suffers from the well-known problem of Gribov copies. Although 
the ambiguity originally noticed by Gribov ||] is not present on the lattice, the maximization procedure used for gauge 
fixing does not uniquely fix the gauge. There are, in general, many local maxima for the algorithm to choose from, 
Oh! each one corresponding to a Gribov copy, and no local algorithm can choose the global maximum from among them. 
<D ■ While various remedies have been proposed pOl HI, they are either unsatisfactory or computationally very intensive. 
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For a recent discussion of the Gribov problem in lattice gauge theory, see Rcf . [[12| . 

An alternative approach is to operate in the so-called Laplacian gauge Jfj|. This gauge is "Landau like" in that 
it has similar smoothness and Lorentz invariance properties |l4]], but it involves a non-local gauge fixing procedure 
that avoids lattice Gribov copies. Laplacian gauge fixing also has the virtue of being rather faster than Landau gauge 
fixing on the lattice. The gluon propagator has already been studied in Laplacian gauge in Refs. (T|, and the 
improved staggered quark propagator in Laplacian gauge in Ref . [fl7| . 

In this report we explore three implementations of the Laplacian gauge and their application to the gluon propagator 
on coarse, large lattices, using an improved action as has been done for Landau gauge in Ref. [||. We study the gluon 
propagator in quenched QCD (pure SU(3) Yang-Mills), using an 0(a 2 ) mean-field improved gauge action. To assess 
the effects of finite lattice spacing, we calculate the propagator on a set of lattices from I0 3 x 20 at f3 = 3.92 having 
a = 0.353 fm to 16 3 x 32 at (3 — 4.60 having a — 0.125 fm. To assist us in observing possible finite volume effects, we 
add to this set a 16 3 x 32 lattice at (3 — 3.92 with a — 0.353, which has the very large physical size of 5.65 3 x 11.30 fm 4 . 

The infrared behavior of the Laplacian gauge gluon propagator is found to be qualitatively similar to that seen in 
Landau gauge. Like Refs. [TJ, [lr| , little volume dependence is seen in the propagator, but, unlike Landau gauge, the 
Laplacian gauge displays strong sensitivity to lattice spacing, making large volume simulations difficult. We conclude 
that further work involving an improved gauge fixing is desired. 



II. THE LAPLACIAN GAUGES 



Laplacian gauge is a nonlinear gauge fixing that respects rotational invariance, has been seen to be smooth, yet 
is free of Gribov ambiguity. It reduces to Landau gauge in the asymptotic limit, yet is computationally cheaper 
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than Landau gauge. There is, however, more than one way of obtaining such a gauge fixing in SU(N). The three 
implementations of Laplacian gauge fixing discussed are 

1. <9 2 (I) gauge (QR decomposition), used by Alexandrou et al. in Ref. |]l5| . 

2. <9 2 (II) gauge (Maximum trace), where the Laplacian gauge transformation is projected onto ST/ (3) by maximiz- 
ing its trace. Also used in Ref. p7[ . 

3. <9 2 (III) gauge (Polar decomposition), the original prescription described in Ref. Jl3| and tested in Ref. fl4|| . 

All three versions reduce to the same gauge in SU(2). 

The gauge transformations employed in Laplacian gauge fixing are constructed from the lowest eigenvectors of the 
covariant lattice Laplacian 

£2 A(C0(z, = X s v(x)l, (1) 

v j 

where 

A(U)(x, = £ [ 2S ( X - y) 5 * j - U^x^Six + A - V) - U,(y)^5(y + p, - x)} , (2) 

i,j = l,...,N for SU(N) and s labels the eigenvalues and eigenvectors. Under gauge transformations of the gauge 
field, 

U^x) -> U^(x) = G(x)U„(x)Gi(x + n), (3) 

the eigenvectors of the covariant Laplacian transform as 

v(x) s -> G(x)v(x) s (4) 

and this property enables us to construct a gauge fixing that is independent of our starting place in the orbit of gauge 
equivalent configurations. 

The three implementations discussed differ in the way that the gauge transformation is constructed from the above 
eigenvectors. In all cases the resulting gauge should be unambiguous so long as the Nth and (N+l)th eigenvectors 
are not degenerate and the eigenvectors can be satisfactorily projected onto SU{N). A complex 2x2 matrix can be 
uniquely projected onto SU(2), but this is not the case for SU(N). Here we can think of the projection method as 
defining its own, unambiguous, Laplacian gauge fixing. 

In the original formulation (l3) , which we shall rather perversely refer to as d 2 (III) , the lowest N eigenvectors are 
required to gauge fix an SU {N) gauge theory. These form the columns of a complex N x N matrix, 

M(x) ij = v(x)i (5) 

which is then projected onto SU(N) by polar decomposition. Specifically, it is possible to express M in terms 
of a unitary matrix and a positive hermitian matrix: M — UP. This decomposition is unique if P = (M^ M) 1 / 2 is 
invertible, which will be true if M is non-singular, i.e., if the eigenvectors used to construct M are linearly independent. 
The gauge transformation G(x) is then obtained by factoring out the determinant of the unitary matrix 

GHx) = U(x)/det[U(x)]. (6) 

The gauge transformation G(x) obtained in this way is used to transform the gauge field (i.e., the links) to give 
the Laplacian gauge-fixed gauge field. G(x) can be uniquely defined by this presciption except on a set of gauge 
orbits with measure zero (with respect to the the gauge- field functional intregral). Note that if we perform a random 
gauge transformation G r (x) on the initial gauge field used to define our Laplacian operator, then we will have 
v(x) s -> v'(x) s = G r (x)v(x) s and M(x) -> M'(x) = G r (x)M(x). We see that P = (M^M) 1 / 2 -> P' = P and hence 
G(x) — > G'(x) = G(x)G\(x). When this is applied to the transformed gauge field it will be taken to exactly the 
same point on the gauge orbit as the original gauge field went to when gauge fixed. Thus all points on the gauge 
orbit will be mapped to the same single point on the gauge orbit after Laplacian gauge fixing and so it is a complete 
(i.e., Gribov-copy free) gauge fixing. This method was investigated for U{1) and SU(2) |14j. It is clear that any 
prescription for projecting M onto some G^(x), which preserves the property G(x) — » G'(x) = G(x)G\(x) under 
an arbitrary gauge transformation G r (x), will also be a Gribov-copy free Laplacian gauge fixing. Every different 
projection method with this property is an equally valid but distinct form of Laplacian gauge fixing. 
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The next approach was used in Ref. and we shall refer to it as <9 2 (I) gauge. There it was noted that only N — 1 
eigenvectors are actually required. To be concrete, we discuss N = 3. First, apply a gauge transformation, G{x) 1 , to 
the first eigenvector such that 

[G(x) 1 v{x) 1 ] 1 = \\v(xY\\ (7) 

and 

[G(x) 1 v(*) 1 ] 2 =[G(x) 1 v(x) 1 ] 3 = V, (8) 

where subscripts label the vector elements, i.e., the eigenvector - with dimension 3 - is rotated so that its magnitude 
is entirely in its first element. Another gauge transformation, G(x) 2 , rotates the second eigenvector, v(x) 2 , such that 

[G{xfv{x)\ = y]{vlf + {vlY, (9) 

and 

[G(x) 2 v(x) 2 } 3 = 0. (10) 

This second rotation is an SU(2) subgroup, which does not act on v 2 (x). The gauge fixing transformation is then 
G{x) = G(x) G(x) . Compare this to QR decomposition, where a matrix, A, is rewritten as the product of an orthog- 
onal matrix and an upper triangular matrix. The gauge transformations are thus like Householder transformations. 

In addition, wc explore a third version, <9 2 (II)gauge, where G(x) is obtained by projecting M(x) onto SU(N) by 
trace maximization. M(x) is again composed of the N lowest eigenvectors and the trace of G(x)M(x) ji is maximized 
by iteration over Cabbibo-Marinari subgroups. 

III. THE GLUON PROPAGATOR IN LAPLACIAN GAUGE 

We extract the gluon field from the lattice links by 

which differs from the continuum field by terms of 0(a 2 ). is then transformed into momentum space, 

M<i) = - E e- lHx+ * /2) Mx + A/2), (12) 

X 

where the available momenta, q, are given by 



A _ 27m M ( _ht hi 



(13) 



is the number of lattice sites in the [i direction. The momentum space gluon propagator is then 

Df v {q) = {Al{-q)A b v {q)). (14) 

Note that this definition includes a factor of Uq 2 from Eq. (]TT|) ; this is the same normalization that was used in 
Ref. §. 

In the continuum, the gluon propagator has the tensor structure plj ] 

DfM - (V - q -^r)6 ab D(q 2 ) + ^6 ab F(q 2 ). (15) 

In Landau gauge the longitudinal part will be zero for all q 2 , but this will not be the case in Laplacian gauge. We 
note that 

q^Df v {q) qu _ xab „ f . 



q 2 



6 ab F(q 2 ) (16) 
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Dimensions f3 a (fm) Volume (fm 4 ) Configurations 



1" 12 3 X 24 460 0.125 1.50 3 X 3.00 100 

2i 16 3 x 32 4.60 0.125 2.00 3 x 4.00 100 

2w 16 3 x 32 5.85 0.130 2.08 3 x 4.16 80 

3 16 3 x 32 4.38 0.166 2.64 3 x 5.28 100 

4 12 3 x 24 4.10 0.270 3.24 3 x 6.48 100 

5 10 3 x 20 3.92 0.353 3.53 3 x 7.06 100 

6 16 3 x 32 3.92 0.353 5.65 3 x 11.30 89 



TABLE I: Details of the lattices used to calculate the gluon propagator. Lattice 2w was generated with the Wilson gauge 
action. 

and use this to project out the longitudinal part. This, unfortunately, makes it impossible to directly measure the 
scalar propagator at zero four-momentum, -D(O). We are, however, able to measure the full propagator, T>(fi). In a 
covariant gauge F(q 2 ) — > |§ where £o is the bare gauge fixing parameter. 

On the lattice the bare propagator is measured, which is related to the renormalized propagator by, 

D(q 2 ) = Z 3 ( fJ ,)D(q 2 ; f x 2 ) F(q 2 ) = Z 3 ( f i)F(q 2 ; ^ 2 ) (17) 

where fi is the renormalization point. In a renormalizable theory such as QCD, the renormalized quantities become 
independent of the regularization parameter in the limit where it is removed. Z 3 is then defined by some renormal- 
ization prescription, such as the off-shell subtraction (MOM) scheme, where the renormalized propagator is required 
to satisfy, 

Dip 2 ) = 4- ( 18 ) 
M 

It follows that 

q 2 D(q)\ q ^=Z 3 (^a). (19) 

With covariant gauge fixing, the longitudinal part is usually treated by absorbing the renormalization into the gauge 
parameter, £o — We shall not discuss the renormalized propagator in this paper, but only consider relative 

normalizations for the purpose of comparing different data sets. 

The ensembles studied are listed in Table 1. To help us explore lattice artifacts, some of the following figures will 
distinguish data on the basis of its momentum components. Data points that come from momenta lying entirely 
along a spatial Cartesian direction are indicated with a square while points from momenta entirely in the temporal 
direction are marked with a triangle. As the time direction is longer than the spatial directions any difference between 
squares and triangles may indicate that the propagator is affected by the finite volume of the lattice. Data points from 
momenta entirely on the four-diagonal are marked with a diamond. A systematic separation of data points taken on 
the diagonal from those in other directions indicates a violation of rotational symmetry. 

In the continuum, the scalar function is rotationally invariant. Although the hypercubic lattice breaks 0(4) invari- 
ance, it does preserve the subgroup of discrete rotations Z(4). In our case, this symmetry is reduced to Z(3) as one 
dimension will be twice as long as the other three in each of the cases studied. We exploit this discrete rotational 
symmetry to improve statistics through Z(3) averaging [|| |). 

As has become standard practice in lattice gluon propagator studies, we select our lattice momentum in accordance 
with the tree-level behavior of the action. With this improved action, 

«-iM¥)+s*'(¥)' (20 » 

this is discussed in more detail in Ref. M. 



IV. RESULTS 
A. Finer lattices 

We start by checking that our finest lattice, at a = 0.125 fm, is "fine enough", by comparing the propagator 
with that of Alexandrou et al. |L5|. Fig. [I] shows the momentum-enhanced propagator, q 2 D(q 2 ), in <9 2 (I) gauge for 
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ensemble 2i (improved action) compared to the data from Ref. jl5| for the Wilson action at ft = 6.0^2j. The two are 
in excellent agreement. 




FIG. 1: The momentum-enhanced propagator in d 2 (l) gauge for the Wilson gauge action at {3 = 6.0 (open circles; data from 
Ref. [jh|) and the improved gauge action at (3 = 4.60 (filled triangles). Both lattices are 16 3 x 32. Data has been cylinder cut. 
The two are in excellent agreement. In Ref. || the difference in normalization between the Wilson and improved actions was 
seen to be ~ 1.08. After taking this into account, the relative normalization here is Z^{j3 — 4.60) = 1.07^3 (/3 = 6.0) 



As the gluon propagator has been extensively studied in Landau gauge, it makes sense to understand the Laplacian 
gauge propagator by comparing it to that in Landau gauge. In accordance with custom, we will discuss the momentum- 
enhanced propagator, q 2 D(q). We define the relative Z% renormalization constant 

Z 3 (Landau) 

Zr = ~^~n — i — : — 7> y 11 ) 
Z3 (Laplacian) 

and choose to perform this matching at fj, = 4.0 GeV. The purpose of this is simply to make the (bare) gluon 
propagators agree in the ultraviolet for easy comparison between the gauges. 

We show the gluon propagator in both Landau and d 2 (II) gauges in Fig. |[ In this figure, the data is for the largest 
finer lattices, 2i and 3. The data has been cylinder cut || || to make comparison easier. As was seen in Ref. j[5| 
the gluon propagator in Laplacian gauge looks very similar to the Landau gauge case matching up in the ultraviolet 
but with a somewhat lower infrared hump. 

Having compared the gluon propagator in 9 2 (II) gauge to Landau gauge, we now compare it to other implementa- 
tions of the Laplacian gauge. We expect each implementation to provide a well defined, unambiguous, but different 
gauge. As we saw when comparing Landau and <9 2 (II) gauges, there is some difference in normalization between the 
propagators in the different gauges, so we define, again at \i — 4.0 GeV 

Z,jd\ll)) 

Zd - w (22) 

In Fig. ||, the momentum-enhanced propagator is plotted in <9 2 (I) and t? 2 (II) gauges for one of the fine lattices 
(ensemble 2i). There is a small relative normalization (Zg — 0.98), but otherwise there is no significant difference 
between them, neither in the ultraviolet nor the infrared. 9 2 (I) and c? 2 (II) also show comparable performance in terms 
of rotational symmetry. 

One difference between Landau and Laplacian gauge is that in the former, the gluon propagator has no longitudinal 
component. We see in Fig. ^ that the longitudinal part of the propagator does indeed vanish in the ultraviolet, 
which is consistent with approaching Landau gauge, but gains strength in the infrared. Comparing 9 2 (I) and c? 2 (II) 
gauges we note that while the transverse parts look alike, the longitudinal behavior is quite distinct. The separation 
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FIG. 2: Comparison of the gluon propagator in Landau and c? 2 (II) gauges for lattice 2i (top) and 3 (bottom) (16 3 x 32, improved 
action, (3 = 4.60 and 4.38 respectively). Data has been cylinder cut. The relative normalizations are Zr{/3 = 4.60) = 1.075 and 
Z R ((3 = 4.38) = 1.20. 



of squares and triangles in c? 2 (II) gauge suggests that F(q 2 ) has stronger volume dependence in that gauge. For a 
comparison between Landau, <9 2 (I) and <9 2 (II) gauges for the quark propagator see Ref. |^7| . 

c? 2 (Ill) gauge works badly, failing even to reproduce the correct asymptotic behavior. Fig. shows data from only 
76 configurations as the gauge fixing failed entirely on four of them. In SU(3), the polar decomposition involves 
calculating determinants which, to our numerical precision, can become vanishingly small, in some cases even turning 
negative at some sites. <9 2 (ILT) gauge was also seen to be a very poor gauge for studying the quark propagator JL7J. 

B. Coarser lattices 

When comparing results from ensembles with different simulation parameters we need to consider three possible 
effects: 
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FIG. 3: The momentum-enhanced propagator from ensemble 2i (/3 = 4.60, 16 3 x 32, improved action) in <9 2 (I) and c? 2 (II) gauges. 
Data has been cylinder cut. The two gauges produce nearly identical gluon propagators, up to a small relative normalization 
(Z g = 0.98). 

1. The dependence of Z% on the lattice spacing, 

2. Errors due to the finite lattice spacing, especially when probing momenta near the cutoff, 

3. Finite volume effects, especially in the infrared. 

In Landau gauge, the dependence of the gluon propagator renormalization, Z^(a), on the cutoff is very weak. Z3 is 
approximately constant with respect to the lattice spacing || . In this case it is easy to compare propagators produced 
on a wide range of lattice spacings. In Fig. ^| we plot the momentum-enhanced propagator on four lattices, which 
have a = 0.125, 0.166, 0.270 and 0.353 fm. We see a very different situation to the one observed in Landau gauge |l: 
the propagators appear to agree in the the deep infrared, yet diverge as the momentum increases. The difference is 
small for the two finest lattices - and non-existent in Fig. |l]- but quite dramatic for the coarsest. 

The correct way to compare the propagators is to normalize them at some common, "safe" momentum, i.e., one 
where we expect finite lattice spacing and finite volume effects to both be small. We choose fi = 0.6 GeV and show 
the results in Fig. [?]. Multiplying the propagator by q 2 in constructing the momentum-enhanced propagator q 2 D(q 2 ) 
reveals a rapid divergence in the ultraviolet. Yet normalizing at higher momenta results in data sets that agree 
nowhere except for at the scale /i. It is interesting that the propagators from ensembles 5 and 6, which have the same 
lattice spacing, have slightly different normalizations. 

Also, as the lattice spacing is increased, the relative normalization between the gluon propagators in d 2 (T) and 
<9 2 (II) gauges slowly diverges from one. As was seen above (Fig. ||), these two Laplacian gauges produce gluon 
propagators with rather different longitudinal components. The longitudinal part of the gluon propagator, multiplied 
by q 2 , is plotted in d 2 (l) and <9 2 (II) gauges for ensembles 4-6 in Fig. |[ using the same normalization determined 
for the transverse parts. Interestingly, the longitudinal part appears to be more affected by the finite volume of the 
lattice than the transverse part. In c? 2 (II) gauge q 2 F(q 2 ) may return to zero as q 2 — > 0, while in <9 2 (I) gauge a small, 
non-zero value appears likely, however, more work is required. 



In previous studies [ fl5[ |l6f it was observed that the infrared gluon propagator saturates at a small volume (~ 1 fm 4 ). 
To further explore this, we also study the propagator at zero four-momentum. As previously discussed, only the full 
propagator, T>(0), can be calculated at zero four-momentum. In order to compare results on all of our lattices we 
normalize the data at 1 GeV. This represents a compromise and is certainly not ideal for all the data sets, hence there 
is some systematic error. These results are shown in Table ||. 

If we restrict ourselves to ensembles 1-5, we see that, given the uncertainties discussed, the sensitivity of V(0) to 
volume is indeed small. We also see the trend, already noted above, for d 2 (l) and d 2 (ll) gauges to become more 
different as the lattice spacing increases, another example of the discretization sensitivity of the Laplacian gauge. In 
the case of the very large lattice, ensemble 6, this sensitivity has become extreme. 
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FIG. 4: Comparison of the longitudinal part of the gluon propagator in d 2 (l) (above) and c? 2 (II) (below) gauges for lattice 
2i (f3 = 4.60, 16 3 x 32, improved action). There have been no data cuts or renormalization. Note that the vertical scales are 
different, so the longitudinal component is smaller in c? 2 (II) gauge than d 2 (l). The data has been sorted into points where all 
the momentum is in the temporal direction, spatial cartestian direction, four-diagonal, and the rest. The separation of temporal 
points (triangles) from spatial cartesian points (squares) suggests that this part of the gluon propagator has stronger volume 
dependence in <9 2 (II) gauge than in d 2 (l) gauge. In a standard covariant gauge this would be a constant, £ (zero, for Landau 
gauge). 

We examine this another way through the transverse propagator, shown in Fig. ^| for ensembles 2i, 3-6. The 
data here corresponds to Table 0, having been normalized at 1 GeV. The propagators are consistent down to low 
momenta, ~ 200 MeV, where we strike trouble. Not only is there a spread between the data sets, but for ensemble 
6, the point from purely temporal momentum is much lower than that from spatial momentum. The situation is the 
same in <9 2 (I) gauge. 





9 



Dimensions p a (fm) Volume (fm 4 ) P(0) - d 2 (l) V(0) - d 2 (ll) 



1 12 3 x 24 460 0.125 lOl 16.6(5) 16.3(4) 

2i 16 3 x 32 4.60 0.125 32.0 17.1(5) 16.0(4) 

3 16 3 x 32 4.38 0.166 99.5 16.7(6) 14.1(3) 

4 12 3 x 24 4.10 0.270 220 19.0(8) 14.6(2) 

5 10 3 x 20 3.92 0.353 300 20.8(9) 14.6(4) 

6 16 3 x 32 3.92 0.353 2040 52(5) 35(2) 



TABLE II: The full propagator at zero four-momentum in d 2 (I) and c? 2 (II) gauges, normalized at 1 GeV for comparison. 



V. CONCLUSIONS 



We have made a detailed study of the gluon propagator on coarse lattices with an improved action in Laplacian 
gauges. We have described and tested three implementations of Laplacian gauge. d 2 (T) (QR decomposition) and 
d 2 (II) (Maximum trace) gauges produce similar results for the scalar transverse gluon propagator, but rather different 
longitudinal components. <9 2 (III) for numerical reasons, works very poorly in SU(3). 

At sufficiently small lattice spacing, the transverse part of the gluon propagator is very similar in Laplacian gauge 
to that in Landau gauge. Laplacian gauge, however, exhibits great sensitivity to the lattice spacing, making results 
gained from coarse lattices difficult to compare with those from finer lattices. This is very different to the situation in 
Landau gauge. By comparing the coarse data sets at sufficiently low momentum, however, we have seen a great deal 
of consistency. In the deep infrared, the results from the largest lattice are difficult to reconcile with the other data. 
Excluding that lattice, the total propagator shows little sign of volume dependence. The most likely explanation of 
the (unimproved) Laplacian gauge results seen here is that on our coarsest lattices (improved lattices with (3 = 3.92, 
4.10, and to some extent even (3 — 4.38), the lattice artifacts are much more severe than in the improved Landau gauge 
case. On the very coarse f3 = 3.92 and 4.10 lattices, it seems very likely that finite volume and discretization errors are 
actually being coupled together by the Laplacian gauge fixing. By implementing an improved Laplacian gauge fixing 
on these lattices we anticipate that these errors will decouple on these lattices and we will be in a better position to 
estimate the infinite volume and continuum limits of the different implementations of Laplacian gauge-fixing. Further 
studies, including an improved Laplacian gauge fixing, will hopefully clarify these issues. 
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FIG. 5: momentum-enhanced propagator from 76 configurations in c? 2 (Ill) gauge from lattice lw ((3 = 5.85, 16 3 x 32, Wilson 
action). There have been no data cuts. This is clearly a very bad gauge fixing. 



7 
6 
5 





= 4.60 
= 4.36 


• 


= 4.10 
= 3.92 




3 

q (GeV) 



FIG. 6: The momentum-enhanced propagator in d 2 (ll) gauge at j3 = 4.60, 4.38, 4.10 and 3.92 (ensembles 2i, 3, 4 & 5). Data 
has been cylinder cut. This figure shows the large sensitivity of Z3 to lattice spacing in Laplacian gauge, quite unlike Landau 
gauge. 



12 



_4 






= 4.10 


/? 


= 3.92 




= 3.92 big 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 

q (GeV) 



3.5 



FIG. 7: The momentum-enhanced propagator from ensembles 3, 4 and 5 in c? 2 (II) gauge. Data has been cylinder cut and 
normalized at 0.6 GeV. 




FIG. 8: Comparison of the longitudinal part of the gluon propagator in <9 2 (I) (top) and d 2 (ll) (bottom) gauges for ensembles 
4, 5 & 6. Data has been cylinder cut. The small and large (3 = 3.92 lattices give diverging results even at large momenta in 
the <9 2 (II) gauge. 
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FIG. 9: Comparison of the transverse part of the gluon propagator in c? 2 (II) gauge for ensembles li, 2-6. Data has been cylinder 
cut. The propagators have been normalized at 1 GeV. 



